Fault-tolerant quantum computation against biased noise 



Panos Aliferis 1 and John Preskill 2 
1 IBM T. J. Watson Research Center, P. O. Box 218, Yorktown Heights, NY 10598 
2 Institute for Quantum Information, California Institute of Technology, Pasadena, CA 91125 

We formulate a scheme for fault-tolerant quantum computation that works effectively against 
highly biased noise, where dephasing is far stronger than all other types of noise. In our scheme, the 
fundamental operations performed by the quantum computer are single-qubit preparations, single- 
qubit measurements, and conditional-phase (cphase) gates, where the noise in the CPHASE gates 
is biased. We show that the accuracy threshold for quantum computation can be improved by 
exploiting this noise asymmetry; e.g., if dephasing dominates all other types of noise in the CPHASE 
gates by four orders of magnitude, we find a rigorous lower bound on the accuracy threshold higher 
by a factor of five than for the case of unbiased noise. 
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Our confidence that large-scale quantum computers 
can be realized has been boosted by the theory of fault- 
tolerant quantum computation [T] , which establishes that 
noisy quantum computers can operate reliably if the noise 
is not too strong. In a fault-tolerant simulation of a quan- 
tum circuit, logical qubits processed by the computer are 
protected from damage using a quantum code, and en- 
coded operations are realized by gadgets constructed from 
the computer's fundamental operations; aside from per- 
forming the desired transformation on the encoded quan- 
tum information, gadgets also exploit the redundancy of 
the code to correct errors caused by the noise. 

Most work on fault-tolerant quantum computation has 
focused on the design of gadgets that work effectively for 
generic noise without any special structure. But actu- 
ally, in many physical settings the noise is expected to 
be highly biased. If the computational basis {|0),|1)} 
coincides with the energy-eigenstate basis for the unper- 
turbed qubit, then typically dephasing (loss of phase co- 
herence in the computational basis, due to entanglement 
with the environment) is far stronger than relaxation (bit 
flips in the computational basis, due to energy exchange 
with the environment). While dephasing arises from low- 
frequency noise, relaxation is dominated by noise whose 
frequency is comparable to the energy splitting. Typi- 
cally, this higher-frequency noise has a different physical 
origin than the low-frequency noise responsible for de- 
phasing, and it can be orders of magnitude weaker. In 
this paper, we analyze fault-tolerant gadgets that are de- 
signed to exploit this bias. 

The fault-tolerant scheme we propose is founded on the 
assumption that the quantum-computing hardware can 
execute a conditional- phase (cphase) gate with highly 
biased noise, where CPHASE is the diagonal two-qubit 
gate with eigenvalues (1,1,1,-1) in the computational 
basis. The complete set of fundamental operations per- 
formed by our quantum computer is 

ftund = {CPHASE, P\ +) ,M a J U {V l+i) ,V lT) } , (1) 

where denotes the measurement of the Pauli op- 

erator cr x , and V\^) denotes the preparation of a sin- 
gle qubit in the state \ip). To construct fault-tolerant 



Gcss operations (see below), we will need to prepare the 
state |+) = (|0) + |1)), and for fault-tolerant universal 
quantum computation, we will also need to prepare the 
states \+i) = ^(| )+t|l)) and \T) = ^(|0) + e«*/*|l)). 
We have not listed the identify operation, which is im- 
plicitly applied whenever a qubit is idle. 

Our central assumption, that the noise in CPHASE gates 
is dominated by dephasing, may apply to some pro- 
posed gate implementations using semiconductor spins 
[5] and superconducting circuits [3] ; it may also apply to 
trapped-ion qubits if the CPHASE gates are driven by mi- 
crowave fields rather than laser pulses [J. Furthermore, 
noise in the preparation of the state |+) is trivially "bi- 
ased" because |+) is an eigenstate of <7 X , and noise in 
the destructive measurement of <r x has no specific struc- 
ture because the measurement has a classical output. We 
show that, through appropriate gadget design, this noise 
bias can be exploited to improve the accuracy threshold 
for quantum computation; e.g., assuming that dephasing 
dominates all other types of noise by four orders of mag- 
nitude, we find that the provable accuracy threshold is 
higher by more than a factor of five than for the case of 
unbiased noise. 

Other authors [5j |6j |7] have proposed fault-tolerant 
gadgets for biased noise, but these previous construc- 
tions work only if the noise is dominated by dephasing 
even for some gates that do not commute with cr z , such 
as the controlled-NOT (cnot) gate. In our view, this as- 
sumption is not physically well motivated. A biased noise 
model should be applicable if, during the execution of a 
gate, the perturbation responsible for the noise couples 
predominantly to the a z components of the qubits. But 
during the execution of a gate that does not commute 
with cr z , the perturbation may not have this property — 
e.g., if the gate is a single-qubit rotation about the x 
axis, then to take into account a possible over-rotation 
or under-rotation of the qubit we should include a per- 
turbation proportional to a x rather than a z . And even if 
the perturbation is dominated by a cr z term, an insertion 
of the perturbation during the execution of a rotation 
about the x axis will propagate to a linear combination 
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of <t z and <7 y , which cannot be described as dephasing 
noise alone. 

Although the biased noise model has a natural formula- 
tion in terms of a Hamiltonian that couples the computer 
to its environment, we will study a stochastic version of 
the model. A stochastic noise model assigns a probabil- 
ity to each fault path — i.e., to each possible set of faulty 
fundamental operations in the circuit. We speak of lo- 
cal stochastic noise with strength e if, for any r specified 
fundamental operations in the circuit, the sum of the 
probabilities of all fault paths with faults at those r lo- 
cations is no larger than e r [8]. In this model no further 
restrictions are imposed on the noise and, in particular, 
the trace-preserving quantum operation applied at the 
faulty locations of each fault path is arbitrary and can 
be chosen adversarially. Therefore, although e quantifies 
the strength of the noise, the faults can be correlated 
both temporally and spatially. It was shown recently in 
[S] that an ideal quantum circuit can be simulated accu- 
rately and with reasonable overhead provided that e is 
smaller than e t h > -67 x 1CP 3 ; this rigorous lower bound 
on the threshold is the best established so far for this 
noise model. 

The noise model that we will analyze in this paper is a 
refinement of local stochastic noise that admits two dif- 
ferent values of the noise strength: e, quantifying the rate 
for faults in preparations and measurements and dephas- 
ing faults in CPHASE gates, and e' <§C e, quantifying the 
rate for all other types of faults in CPHASE gates. In this 
model, a fault path indicates not only which locations 
are faulty, but also, for each CPHASE gate, whether a de- 
phasing fault or some other type of fault has occurred. 
We speak of local stochastic biased noise if the sum of the 
probabilities of all fault paths that are faulty at r spec- 
ified locations, of which s are non-dephasing faults at 
CPHASE gates, is no larger than e r ~ s (e') s . For dephasing 
faults, all Kraus operators are assumed to be diagonal in 
the computational basis, and for all other types of faults, 
the Kraus operators are arbitrary. We refer to the ratio 
e/e' as the noise "bias." 

Our scheme for fault-tolerant quantum computation 
will be protected by a code C\ > C2 where > denotes code 
concatenation. The inner code C\ protects against de- 
phasing, mapping highly biased noise to a more balanced 
effective noise model with reduced noise strength. The 
code C2 maps unstructured noise with strength below eth 
to noise with negligible strength. We take C\ to be a 
lcngth-n repetition code in the dual basis, where n is 
odd; the n—1 check operators are J®- 7 <g)cr x <g)er x <g)/' 8m_: '~ 2 
(j = 0, 1, . . . , n — 2), and the logical Pauli operators act- 
ing on the one encoded qubit are er x = <7^®I®I® - ■ -(g)/, 
and a% — cr z ® cr z <£> • • • <£> er z . The code C\ can correct 
(n — l)/2 er z errors but provides no protection against 
a x errors. We take C2 to be a concatenated CSS code. 
For a CSS code [8] , the fault-tolerant encoded versions of 
operations in the set 

Gcss = {CNOT, P\ 0) , V\ +) , M a . , M a J (2) 
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FIG. 1: Scheme for achieving fault-tolerant universal quan- 
tum computation against highly biased noise. 

can be built out of operations that are contained in this 
set; furthermore, Gcss operations suffice for measuring 
the error syndrome. 

We will use the fundamental operations in Giund (where 
the only necessary state preparation is to construct 

Gcss gadgets protected by C%. Combining with known 
constructions for CSS codes [TU], we obtain Gcss gad- 
gets protected by C\ >C2- Finally, CSS operations will be 
extended to a universal set by appending preparations 
of the states \+i) and \T); high fidelity encoded copies 
of these states can be prepared by teleporting ( "inject- 
ing" ) into the C\ > C2 block and then performing state 
distillation [11 . Our scheme for achieving fault-tolerant 
universal quantum computation is illustrated in Fig. [I] 
(we denote an operation O encoded in C\ or C\ > C2 as 
O l or O respectively). 

The crux of our gadget constructions, and the basis for 
our threshold analysis, is the design of the Ci-protected 
CNOT gadget using the operations in Gtund- The key idea 
is to use a variant of state teleportation that simultane- 
ously executes the encoded gate and extracts the error 
syndrome. But first, let us discuss how to construct C\- 
protected gadgets for the other operations in Gcss- 

The destructive measurement of cr x is performed by 
measuring er x for all qubits in the code block, and then 
computing the majority of the measurement outcomes. 
This measurement is fault tolerant in the following sense: 
if m of the qubits in the input code block have errors, and 
s of the single-qubit measurements are faulty, then the 
outcome of the noisy encoded measurement agrees with 
the outcome of an ideal encoded measurement provided 
that m+s < {n— 1)/2. The preparation of is exe- 
cuted transversally: V\+) L = {^1+))^ ■ This operation 
is fault tolerant because at least (n+l)/2 of the prepara- 
tions of |+) must be faulty to cause an encoded error. 

A nondestructive measurement of is performed with 
the circuit depicted in Fig. [2j an ancilla qubit is pre- 
pared in the state |+), n consecutive CPHASE gates are 
executed, and then <r x is measured on the ancilla qubit. 
If performed only once, this measurement is not fault tol- 
erant, because a single cr z error acting on the ancilla can 
flip the outcome; however, fault tolerance can be ensured 
by repeating the measurement r times, where r is odd, 
and computing the majority of the outcomes. Although 
<r z errors acting on the data qubits do not affect the mea- 
surement outcome, they might contribute to a logical 
error that could affect subsequent operations. It is there- 
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FIG. 2: Gadget that measures (here, n = 3). The mea- 
surement is repeated r times to ensure fault tolerance (here, 
r = 3), with the repetitions staggered as shown so that the 
data qubits are never idle in between consecutive gates. 
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FIG. 3: On the left, a "one-bit teleportation" circuit. On the 
right, the error correction gadget. 



fore noteworthy that if the input block has m cr z errors 
and the measurement gadget has s dephasing faults, then 
there will be no more than m+s <r z errors in the output 
block. 

The preparation of |0)x, is executed by first preparing 
\+)l, and then performing the nondestructive measure- 
ment of a z . Again, <r z errors acting on the data qubits 
will not disturb the eigenvalue of a z \ only the <r z er- 
rors acting on the ancilla are problematic. Therefore, 
the fault tolerance of M. a L ensures the fault tolerance of 
7-'|o) i , ■ If the measurement result is a z = —1, then the 
prepared state differs from |0)l by a known logical 
error. This error need not be corrected; rather it is used 
to update the "Pauli frame" of the encoded block [P2"] . 

Error correction can be performed by teleporting an 
encoded block [12]. Because we only need to correct ct z 
errors, the error correction gadget can be simplified to an 
encoded version of the "one-bit teleportation" circuit [T3] 
depicted in Fig. [3j Ideally, the output encoded qubit has 
the same state as the input encoded qubit, apart from a 
possible a z (if the outcome of M. a L is —1) and a possi- 
ble ct^ (if the outcome of Ai a L a L is —1); thus the Pauli 
frame is updated based on the measurement outcomes. 
The nondestructive measurement of <J z a z is performed 
using one ancilla qubit and 2n CPHASE gates, and is re- 
peated r times, much as for the measurement of a z de- 
scribed above. If the input block has m <r z errors and 
the error-correction gadget has s dephasing faults, then 
the outcome of M a L agrees with the ideal case, provided 
m+s < (n— 1)/2. Furthermore the outcome of M a L a L 
agrees with an ideal measurement for s < (r— 1)/2, and 
the number of cr z errors in the output block is no more 
than s. 
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FIG. 4: Fault-tolerant CNOT gadget. Pauli operators that 
update the Pauli frame (not shown) are determined by the 
measurement outcomes. 



By combining one-bit teleportation gadgets acting on 
both output blocks with a logical CNOT gate, we obtain 
the CNOT gadget depicted in Fig. [3] where the upper 
block is the control block and the lower block is the 
target block. The measurements of <r z a z and a z a z a z 
are repeated r times using ancilla qubits and CPHASE 
gates, as described previously. If the input target block 
has mi a z errors, the input control block has m-i er z er- 
rors, and the CNOT gadget contains s dephasing faults, 
then both M. a L 's agree with the ideal case provided that 
mi+s < (n— 1)/2 and m2+s < (n— 1)/2; furthermore, 
each output block contains no more than s a z errors, 
and M. a L a L and M a L a L <7 L agree with the ideal measure- 
ments for s < (r— 1)/2. Further properties of the CNOT 
gadget are discussed in Appendix A. 

Now consider a circuit constructed from these C\- 
protected Qcss gadgets — how accurately does it simulate 
an ideal circuit? A gadget operates correctly if all of the 
encoded measurements it contains agree with the ideal 
case (the case in which the input blocks have no errors 
and the gadget contains no faults); otherwise the gadget 
fails. For each (Jess gadget, we have derived an upper 
bound on its probability of failure in terms of the noise 
strength e and the bias factor e/e' of the local stochastic 
biased noise model. See Appendix B for details of this 
combinatorial analysis. 

The largest of these upper bounds (found for the CNOT 
gadget) is denoted e^; it can be regarded as the effec- 
tive noise strength of a local stochastic noise model that 
characterizes the noise in Ci-protected Gcss circuits. In 
Fig. [5] we have plotted as a function of e for two dif- 
ferent values of the bias. If is below the previously es- 
tablished lower bound on the accuracy threshold for CSS 
operations (e^f s > .67 x 10 -3 [9]), then we can choose 
the code C 2 so that the Gcss gadgets protected by C\ >C 2 
are arbitrarily accurate. Thus we set e^ 1 ^ = .67 x 1CP 3 , 
and choose n so that e is as large as possible. If the 
bias is 10 4 , then the maximum value is e — 2.50 x 1CP 3 , 
occurring at n = r = 1 1 . 

Furthermore, as we show in Appendix C, the injection 
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FIG. 5: Upper bounds on the effective noise strength e' 1 ' for 
Ci-protected 5css operations as a function of e and the bias 
e/e' (where for each value of e, we optimize over n and r). 
The straight line with slope unity serves as a guide to the eye. 

and distillation of the \+i) and |T) states can be per- 
formed effectively for e < 2.50 x 10 -3 . We conclude that 
2.50 x 10 -3 is a lower bound on the accuracy threshold 
for universal quantum computation under local stochas- 
tic biased noise with bias 10 4 . This is an improvement 
by about a factor of four compared to the case of un- 
structured noise. On the other hand, the improvement 
is more modest for smaller values of bias; e.g., for a bias 
of 10 3 , our lower bound becomes 1.54 x 10~ 3 (where the 
maximum value e — 1.54 x 10~ 3 occurs at n = r = 7), 
an improvement by about a factor of two. 

These results can be further improved by modifying 
the decoding procedure for C\ > C2 ■ Recall that the out- 
come of each Ci-protected measurement is determined 
by a majority vote, and note that this outcome is more 
likely to be wrong if the vote is "close" — i.e., if the ma- 
jority has just one more vote than the minority. For 
example, if (n— 1)/2 of the a x measurements inside M. a L 



disagree with the choice of the majority, then there might 
be (n+l)/2 errors in the block, resulting in an encoded 
error. But if only (n— 3)/2 qubits disagree with the ma- 
jority, then there must be at least (n+3)/2 errors to cause 
an encoded error. This observation also applies to the 
majority voting in A4 (J l !T l and M. rT t a L (7 L. 

The code C2 can be decoded more reliably by exploit- 
ing information concerning which Ci-protected measure- 
ments have close votes; see Appendix D for details. Using 
this more sophisticated decoding method, we find that 
the accuracy threshold improves to 2.09 x 10" 3 for bias 
10 3 , and improves to 3.51 x 10 -3 for bias 10 4 . 

Fault-tolerant methods will be essential for achieving 
large-scale quantum computation. These methods can be 
more effective when customized for the particular prop- 
erties of the noise in the computing hardware. In this pa- 
per, we have explained how to exploit noise asymmetry in 
diagonal gates to make fault-tolerant quantum comput- 
ing work better. We have analyzed the performance of 
our scheme for local stochastic biased noise; using tech- 
niques described in [HI Q31 [15] , a more realistic Hamilto- 
nian biased noise model could also be analyzed. 

A notable property of our constructions is that the only 
fundamental operation used by the quantum computer, 
other than single-qubit preparations and measurements, 
is the two-qubit CPHASE gate. This feature is attractive, 
because in some physical settings CPHASE gates are par- 
ticularly easy to execute with highly biased noise and 
reasonable fidelity; e.g., our companion paper |16j dis- 
cusses how the encoding scheme we have formulated here 
applies to superconducting flux qubits. 
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92-ER40701, NSF under Grant No. PHY-0456720, and 
NSA under ARO Contract No. W911NF-05-1-0294. 
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APPENDIX A: 



COMMENTS ON THE CNOT 
GADGET 



In our discussion in the main text, the measurement of 
cr^a^ is repeated r times inside Ai a L a L and the mea- 
surement of a^a^a^ is also repeated r times inside 
M rT L a L rT L . In general, the number of repetitions could 
be different in the two cases, 7*1 and for the measure- 
ments of a^a^ and o^o^o^ respectively, and for now we 
will distinguish r\ and r2 from n so that the counting we 
describe below will be more transparent. In fact, later 
on we will set r\ — T2 — n, which turns out to be opti- 
mal or nearly optimal in the cases we have studied. A 
further advantage of the choice r\ — — n is that we 
can eliminate storage locations (where quints are idle) in 
the CNOT gadget by staggering the measurements as in 
Fig. [2j for this reason we will not include any faults at 
storage locations in our estimate of the failure probabil- 
ity. 

A noteworthy property of the CNOT gadget is that, if 
ri = T2 = n, and if the measurements are staggered as in 
Fig. [2] then the latest operations on the output blocks act 
one time step before the earliest operations on the input 
blocks. This property is obscured by the diagrammatic 
notation in Fig. [4] but it is evident once we consider the 
full circuit as in Fig. [6] Let us say that a data qubit "in- 
teracts" in a time step in which it is coupled to an ancilla 
qubit by a CPHASE gate. We choose a standard ordering 
for the n qubits in each block, such that the interactions 
of qubit j lag j—1 time steps behind the interactions of 
qubit 1. Then, in the CNOT gadget, qubit 1 in the out- 
put control block interacts during time steps 1 through 
n, and qubit 1 of the output target block interacts during 
time steps 1 through n. Meanwhile, qubit 1 in the input 
control block interacts during time steps n+1 through 
3n, and qubit 1 in the input target block interacts dur- 
ing time steps n+1 through 2n. Therefore, in time step 
n+1, as qubit 1 in each input block begins to interact, 
qubit 1 in each output block is ready for execution of 
the next gadget. This is a characteristic feature of gate 
teleportation; it implies that a circuit of Gcss operations 
can be simulated in constant depth, independent of the 
size and depth of the simulated circuit. 



APPENDIX B: THE THRESHOLD FOR g C ss 
OPERATIONS 

Among all Gcss gadgets, the CNOT gadget contains the 
largest number of fundamental operations. Therefore, if 
we derive an upper bound on its probability of failure, 
this bound will also apply to all the other Gcss gadgets. 

To estimate the probability of failure for a CNOT gad- 
get, we first observe that a tr x error on a data qubit or a z 
errors acting on many ancilla qubits can cause an incor- 



rect outcome of M.„l„l or M. C7 l c .l c 



The outcome of 



M. a L may be incorrect due to <r z errors acting on multi- 
ple qubits in a single block, or due to a cr x error acting on 
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FIG. 6: The full circuit of the CNOT gadget (here, for com- 
pactness we have chosen r\ = r-2 — n = 2, but the structure 
of the circuit is similar for odd values). 



an ancilla qubit that propagates repeatedly to generate 
many er z errors in the code block. 

In addition, we must take into account errors in preced- 
ing gadgets that could propagate into the CNOT gadget 
we are considering; e.g., Fig. [7] depicts a CNOT gadget 
preceded by CNOT gadgets acting on each of its input 
blocks. A cr z error in one of the preceding gadgets might 
affect the outcome of aM^ in the current gadget, and 
a <7 X error in one of the preceding gadgets can alter the 



outcome of A4, 



or M r . 



in the current gadget. 



To understand the effect of tr x errors in the CNOT gad- 
gets, note that an incorrect outcome of M^l^l produces 
a <7 X error acting on the CNOT gadget's output control 
block, and an incorrect outcome of M (y L a L a L produces 
a cr x error acting on the CNOT gadget's output target 
block. Consider, e.g., the control block of the final CNOT 
gadget for the case depicted in Fig. [7] and suppose that 
a single non-dephasing fault in a CPHASE gate contained 
in the AAuL^l^l of the immediately preceding gadget al- 
ters the outcome of that measurement and also of the 
A4 a L a L in the later gadget. Then this one fault causes 
logical errors in each of two consecutive gadgets. How- 
ever, the logical error in the earlier CNOT gadget is a <r x 
acting on its output target block, and has no effect on 
its output control block. Therefore, we can propagate 
the logical error forward from the earlier gadget to the 
later gadget; i.e., we may just as well say that the earlier 
gadget is executed properly, and the logical error occurs 
only in the later gadget. 

More generally, whenever a single er x error causes log- 
ical errors in two consecutive gadgets, the logical error 
in the earlier gadget can be propagated forward into the 
later gadget in this way. Thus, we may hold the o~ x error 
responsible for only the failure of the later gadget, and 



FIG. 7: A CNOT gadget preceded by CNOT gadgets acting on each of its input blocks. 



we may disregard the damage it inflicts on the earlier 
gadget. In effect, then, a single non-dephasing fault oc- 
curring with probability e' can cause the failure of only 
one of our Ci-protected gadgets. 

A measurement of cr^a^ uses 2n CPHASE gates and a 
measurement of a^a^a^ uses 3n CPHASE gates. There- 
fore, the CNOT gadget contains (2n + 3r 2 )n CPHASE 
gates. We pessimistically assume that any non-dephasing 
fault in a CPHASE gate that is either contained within the 
CNOT gadget or that propagates into the CNOT gadget 
causes the gadget to fail. We denote by e^j the proba- 
bility of failure due to a non-dephasing fault in a CPHASE 
gate and conclude that 



< (2r 1 + 3r 2 + 2r)ne' 



(Bl) 



where r < max(ri,r2). Here, for each input block, rn 
is an upper bound on the number of CPHASE gates in 
the preceding gadget where errors can propagate into the 
current gadget. 

Now, suppose that there are no CPHASE gates with 
non-dephasing faults, and consider the probability of 
failure of the CNOT gadget due to dephasing faults at 
CPHASE gates, and due to faults in preparations and mea- 
surements. We may assume without loss of generality 
that the faults at the operations and A4 a ^ have di- 
agonal Kraus operators, as cr x errors have no effect at 
these locations. 

A dephasing fault can alter the outcome of a cr^a^ or 
a^cr^a^ measurement only if it acts on the ancilla qubit 
used during the measurement. For each of these logi- 
cal measurements there is one preparation and one mea- 
surement operation; furthermore, there are 2n CPHASE 
gates for the cr^a^ measurement and 3n CPHASE gates 
for the <Jz (J z (J z measurement. We therefore obtain up- 
per bounds on the probability of failure for the majority 



vote of the repeated measurements: 



! ( r J i j((2n + 2) £ ) n ^ , (B2) 
e{M aiaivi ) < I r r 2 2 +1 ) ((3n + 2) e)^ . (B3) 



A measurement of cr^ can fail if the majority of the 
qubits in the block have errors. For each qubit, the error 
can arise during the preparation of the qubit, the mea- 
surement of the qubit, or a CPHASE gate that acts on the 
qubit. Therefore, an upper bound on the probability of 
failure is 



e{M [ l\) < ((r + r 1 + r 2 + 2)ey 



for the measurement of the control block, and 



e(M l *l)<( n n +1 )((r + r 2 + 2)e)^ 1 



(B4) 



(B5) 



for the measurement of the target block. Denoting by 
e)P the probability of failure due to faults other than 
non-dephasing faults in CPHASE gates, we obtain 



< 2 



n+l 
1 



(271 + 2)^ + (3n + 2)^ £ ,A r (B6) 



where to obtain the second inequality we have substi- 
tuted T\ = r-2. = T = n. Our upper bound on the total 
probability of failure for the CNOT gadget is 



(B7) 
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FIG. 8: Teleportation circuits for the Clifford group gates Q = 



O i(ir/4)o 



and S 

-i(7r/8)cr z 



and for the non-Clifford gate 
T = e -1 ^/ 8 )^ p or t ne Clifford group gates the measurement 
determines a Pauli operator that updates the Pauli frame, 
and for the non-Clifford gate the measurement determines a 
non-Pauli correction that must be applied in the next step. 



The quantity is the effective noise strength for our 
Ci-protected Gcss gadgets. In particular, if we set r = 
n = 11, e = 2.50 x 1CT 3 , and e/e' = 10 4 , we find < 
.67 x 10~ 3 so that the effective noise strength is below the 

-CSS 



threshold eJ^ SM > .67 x 10~ 3 for Gcss gadgets protected 
by C2 [5]. Thus 2.50 x 10~ 3 is a lower bound on the 
accuracy threshold for Gcss operations assuming a local 
stochastic biased noise model with bias 10 4 . 



APPENDIX C: ACCURACY THRESHOLD FOR 
UNIVERSAL QUANTUM COMPUTATION 

If e is below the Gcss threshold, Gcss gadgets pro- 
tected by C\ > C2 are highly reliable. To extend our gad- 
gets to a universal set protected by C\ > C2 , we use gate 
teleportation as shown in Fig. [8] Provided we can pre- 
pare the state \+i), we can use the Gcss operations cnot, 
"P\+s, A^ CTx , and M. Cl to teleport Q = exp(i(7r/4)cr x ) and 
S = exp(— i(7r/4)cr z ). Together with the cnot gate, Q 
and S suffice to generate the Clifford group. Provided we 
can prepare the state \T), we can go beyond the Clifford 
group and achieve universality by using Gcss operations 
and S to teleport the gate T — cxp(— i(ir/8)a z ). Thus, we 
can do reliable universal quantum computation if we can 
perform CSS operations reliably and we can also prepare 
high-fidelity copies of the state \+i) (the eigenstate of a y 
with eigenvalue +1) and the state \T) (the eigenstate of 
Sa x with eigenvalue +1). 

Furthermore, if we are able to prepare noisy copies of 
\+i) and \T) that are not too noisy, then high fidelity 
copies can be generated via state distillation protocols 
These protocols are based on CSS stabilizer codes, 
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FIG. 9: Injection of the state |i/>) into the C\\>Ci block by using 
teleportation. After the encoded Bell state |$o) is prepared, 
one Ci >C2 block is decoded to C\. Then, a Bell measurement 
is performed on the Ci block and the input state \ip). 



for which Gcss operations suffice to measure the error 
syndrome and to decode. The distillation protocol for 
uses Steane's [[7,1,3]] CSS code. In each round of 
the protocol, the code's check operators are measured 
for seven noisy copies of the input state; the encoded 
qubit is accepted and decoded if the error syndrome is 
trivial. The state \+i) is prepared successfully unless at 
least three of the input states have errors. Similarly, the 
distillation protocol for \T) uses a [[15,1,3]] CSS code. In 
each round of the protocol, the code's check operators 
are measured for fifteen noisy copies of the input state; 
the encoded qubit is accepted and decoded if the error 
syndrome is trivial. Here, too, the state |T) is prepared 
successfully unless at least three of the input states have 
errors. 

The error threshold for the |T) distillation protocol was 
estimated in |llj . where it was shown that an input er- 
ror probability as high as 14.1% can be tolerated if each 
input state is "twirled" by applying Sa x with probabil- 
ity 1/2. The error threshold for distillation is even 
higher. Therefore, if Gcss gadgets protected by C\ > C2 
are reliable, and we can also inject input states into the 
C\ \>C2 block with probability of error below 14.1%, then 
we can do reliable universal quantum computation. (By 
distilling \+i), we can teleport the S gate, enabling us to 
perform the twirling step in the |T) distillation protocol.) 

The state injection is performed by teleportation as 
in Fig. [9] Let us use \tp) to denote a state encoded in 
C\ >C2, to distinguish it from the state encoded in 

C\. To inject the single-qubit state \ip) into the C% > C2 
block, first the encoded Bell state |$ ) = -^(|Q)|0) + 
|1)|1)) is prepared, and then one of the code blocks is 
decoded to C\. To complete the teleportation, a joint 
Bell measurement is performed on the C\ block and the 
unprotected state This procedure prepares the state 
I i/j), up to a logical Pauli operator that is known from the 
outcome of the Bell measurement. 

Because |$o) can be prepared by using Gcss gadgets 
which arc well protected by C\ \>C2 (we start with the en- 
coded state |+) and the encoded state |0), and we apply 
an encoded CNOT gate), we may assume that the prepa- 
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ration of |$o) 1S flawless. Thus the state injection might 
fail because of a decoding error, because of an error in 
the Bell measurement, or because of a fault during the 
single-qubit preparation of \tp). 

Now let us suppose that C2 = C >k is obtained by con- 
catenating the CSS code C all together k times. The 
decoding of C\ > C2 is performed recursively: In the first 
step, Ci >C >k is decoded to C\ > C^' -1 ) using gadgets 
protected by d t>C >( - k -^; then, C x > C 1 ^- 1 ) is decoded 
to Ci > C >( - k ~ 2 ) using gadgets protected by C\ > C >( - k ~ 2 \ 
and so on. In the last step, C\ >C is decoded to Ci using 
gadgets protected by C\. Let us denote by e(V) the prob- 
ability that a logical error occurs at any level during this 
recursive decoding. If decoding is staggered so that no 
qubits are idle during the Bell measurement, and if the 
measurement of <r z crf' is repeated r times inside A4 a ^ a L, 
the probability of a state injection error is 



(CI) 



where e accounts for the probability of a fault in the 
single-qubit preparation of and 



< (2rn + r)e' 
+ U0((nH 



"(1- 
3) £ ) 



- r)e 

r+l 



t )((2r + 2) £ ) 5 



Here, (2rn+r)e' bounds the probability of error in the 
Bell measurement due to a non-dephasing fault in a 
CPHASE gate; the fault could occur in one of the r(n+l) 
gates contained in M. a ^ a L, or in one of the rn gates con- 
tained in a measurement in the immediately preceding 
Ci-protected CNOT gadget (which is part of the recursive 
decoding circuit). Furthermore, (l+r)e bounds the prob- 
ability of error in A / l (Tx due to a dephasing fault; the mea- 
sured qubit participates in r CPHASE gates contained in 
M. a ^ a L and also in the measurement itself (the probabil- 
ity of a fault in the preparation of this qubit in the state 
\tp) has already been included in Eq. (jClJ). The next 
to last term bounds the probability of error in M a a L 
due to a dephasing fault; the ancilla qubit in each of the 
r measurements inside M ai<J L participates in one "P\+), 
ri+1 CPHASE gates, and one M. 0x . Finally, the last term 
bounds the probability of an error in M a L ; each qubit in 
the measured block participates in one P\+), r CPHASE 
gates contained in the preceding Ci-protected CNOT gad- 
get, r CPHASE gates contained in M. a a L, and one M ax . 

If we set r = n = 11, e < 2.50 x 10~ 3 , and e' < 1CT 4 £, 
we find e BM < 3.01%. Since e(V) < 8.24% for £« < 
.67 x 10~ 3 jjj], we conclude that £(P|^) < 11.5% which is 
below the 14.1% distillation threshold. Thus 2.50 x 10~ 3 
is a lower bound on the accuracy threshold for universal 
quantum computation under local stochastic biased noise 
model with bias 10 4 . 
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FIG. 10: A Bell state |$o)l is prepared by starting with two 
blocks in the state \+}l and then measuring a^a^. 



APPENDIX D: IMPROVED THRESHOLD VIA 
FLAGGING AND MESSAGE PASSING 



We can improve our lower bound on the accuracy 
threshold by using a more sophisticated decoding proce- 
dure for C\ >C2- We note that the syndrome information 
for C\ is helpful for optimizing the decoding of C2 in the 
concatenated block; yet the decoding procedure that we 
have described so far makes no use of this information — 
after C\ is decoded, the Ci-syndrome is discarded. Now 
we consider a new decoding procedure where some infor- 
mation about the Ci-syndrome is retained and used in 
the decoding of C2 ■ 

For the sake of clarity, we continue to make a distinc- 
tion between ri, T2 and n, even though we will set them 
equal later on. We say that a vote is "close" if the win- 
ners have one more vote than the losers. Thus M. a L has a 
close vote if the error syndrome indicates (n— 1)/2 errors 
in the block, A4 a L a L has a close vote if (r\ — 1)/2 of the 
measurements disagree with the majority, and Ai a L a L a L 
has a close vote if (v2— 1)/2 measurements disagree with 
the majority. If a gadget contains no close votes, then 
we decode C\ as usual. But if the gadget contains a close 
vote, then a flag is raised after decoding. The flag signi- 
fies that the gadget has a higher than usual probability 
of failure, information that will be exploited during de- 
coding at the next level up in the concatenated block, 
using a scheme described in [j"2"] . 

For simplicity, we consider a version of the scheme in 
[12] where C2 is the concatenated 4-qubit code (C4) with 
check operators <7® 4 and <r® 4 ; this is the case analyzed 
in [5]. The basic building blocks for the construction 
of C2-protected Gcss operations in this scheme are Bell 



states |$ )= 75(|0)|0) + |1)|1)) as on the left of Fig. 10 
and CNOT gates followed by single-qubit measurements 
as on the left of Fig. [IT] We can then construct encoded 
versions of these two basic building blocks by using C\- 
protected gadgets as on the right of Figs. [TO] and [TT] 

Consider how to refine our estimate of the failure prob- 
ability for Ai rT L a L in Fig. 10 taking flagging into account. 



If there is a flag, then (ri+l)/2 faulty measurements 
might cause M. a L a L to fail, but without a flag at least 
(ri+3)/2 faulty measurements are required. Therefore, 
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FIG. 11: On the left, the CNOT gates and measurements used to implement the teleportations of C2 blocks and subblocks in 
Fig. 3 in On the right, a Ci-protected implementation of the operations on the left. The outcome of a measurement of 
<Tz(Tz or a^a^a^ on the right corresponds to the outcome of the measurement of a z on the left with the same superscript, 
while the outcomes of the measurements of cr x on the right are combined as shown to give the outcomes of the measurements 
of <T X on the left. 



instead of Eq. ( B2 1 we have 



e f (M„ iai ) < (^j ((27^2) £ )^ 
for the probability of failure with a flag, and 



,f(M a L a L) < 



ri+3 



((2n + 2)e)" 



Cl±3 



(Dl) 



(D2) 



for the probability of failure without a flag. In fact, be- 
cause the circuit in Fig. [TO] prepares an ancilla state, we 
find it advantageous to repeat this measurement of cr^cr^ 
a smaller number of times t < n and to postselect on the 
cases without a flag. In this case, the conditional proba- 
bility that M a L a L fails when it is accepted is 



value is achieved by choosing r± = r% = n = 7 and t = 5). 
In addition, by an analysis similar to the discussion in 
Appendix C, we find that for e = 3.51 x 1CP 3 the prob- 
ability of an error in state injection is eCPu\) < 10.4%, 
which is below the 14.1% distillation threshold. Thus 
3.51 x 1CP 3 is our improved lower bound on the accuracy 
threshold for universal quantum computation under local 
stochastic biased noise with bias 10 4 . 



z{M a L a L\a,cc) < 



(iU) ((2n + 2)e)^ +2nte' 



1 



' t N 

t-1 . 



/ ((2n + 2)e) 2 - 2nte' 

(D3) 

t \ — - 
where (t-i) ((2n + 2)e) 2 bounds the probability that 

2 

a flag is raised, and 2nte' is the probability of a non- 
dephasing error in one of the CPHASE gates. 

We can perform a similar analysis to bound the proba- 
bility of failure, with and without a flag, for the measure- 
ments M. a L a L, M. a L a L. a L, and M.„h in Fig. fTTj These 
bounds can now be plugged into the analysis 01 the de- 
coding of C 2 in M- For bias e/e' = 10 4 , we find that 
reliable {C\ > C2)-protected Gcss operations can be im- 
plemented if e is below 3.51 x 10 -3 (where this optimal 



